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ABSTRACT 

We present measurements of carbon monoxide emission in the central region 
of the nearby starburst NGC 6000 taken with the SubmiUimeter Array. The J — 
2-1 transition of ^^CO, ^^CO, and C^^O were imaged at a resolution of ~ 3" x 2" 
(450 X 300 pc). We accurately determine the dynamical center of NGC 6000 at 
Q;j2ooo.o = 15'^49'"4g:5 and 5j2ooo.o = -29°23'13" which agrees with the peak of 
molecular emission position. The observed CO dynamics could be explained in 
the context of the presence of a bar potential affecting the molecular material, 
likely responsible for the strong nuclear concentration where more than 85% of the 
gas is located. We detect a kinematically detached component of dense molecular 
gas at relatively high velocity which might be fueling the star formation. A 
total nuclear dynamical mass of 7 x 10^ Mq is derived and a total mass of gas 
of 4.6 x IO^Mq, yielding a Mgas/Mdyn ~ 6%, similar to other previously studied 
barred galaxies with central starbursts. We determined the mass of molecular gas 
with the optically thin isotopologue C^^O and we estimate a CO-to-H2 conversion 
factor Xco = 0.4 x 10^^ cm~^(K kms~^)~^ in agreement with that determined in 
other starburst galaxies. 

Subject headings: galaxies: fundamental parameters — galaxies: ISM — galaxies: 
kinematics and dynamics — galaxies: nuclei — galaxies: starburst — radio lines: 
galaxies 
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Introduction 



The relation between tlie morphology and kinematics of the gas in the inner region of 
galaxies and the connection to its nuclear activity is not completely understood. Nuclear 
regions are often obscured by dust at optical wavelengths but can be viewed in continuum 
infrared dust emission as well as molecular line emission. CO can be used as a tracer of 
molecular gas, and its rotational transitions are observed at millimeter and sub-millimeter 
wavelengths. Observations of the morphology and kinematics of the molecular gas content 



in the cent ral regions q : 



there (e.g., 



Jo gee et al 



spirals can provide insight into the dynamical processes occurring 



2005 



Perez-Ramirez et al. 



20001 ). Studies indicate that instabilities 



in bars and dissipation of gas clouds remove angular momentum from material orbiting the 
galactic center, driving gas and dust inward to fuel star formation 



1990 



Sakamoto et al. 



1999 



Knapen et al. 



2002 



Jogee et al. 



2005 



e.g., 


Pfennieer & 


N^orman 


Sheth et al. 


2005 


)• 



NGC6000 is a nearby (D~31.6Mpc, 



Pizzella et al. 



20051 ) barred spiral starburst galaxy. 



Despite its proximity and brightness, this galaxy has not been studied in as great detail as 
others of similar distance and luminosity, due to its southern declination. Hubble Space 



Telescope observations reveal 



a large-scale (> 1 kpc) bar (ICaroUo et al. 



several brigh t sources in circumnuclear star-form ing rings and 



1997 



2002 



Fathi fc Peletier 



2003|). Ho wever, 



no such barred structures a re identified towards the more irregular nuclear region (jCaroUo 



1999 



Hunt fc Malkan 



20041 ). Infrared observations towards the nucl eus of NGC 6000 hav e 



found significant polycyclic aromatic hydrocarbon (PAH^ 



20041 ) and determined a dust ter nperature of 29.2 ±2.6 K (lYang &: Phillips 



appe ars bright in CO emis sion ( lYoung et al. 



1995 



emission (ISiebenmorgen et al. 



H I f iKoribalski et al. 



Mauersberger et al.l 



20071). NGC 6000 



19991) as well as 



20041 ) as detected with single-dish telescopes. However no detailed 
morphological and kinematical study of the gas content has been performed yet. A 
compilation of the observed and derived properties for NGC 6000 are presented in Table [1] 
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In this paper we present observations of the carbon monoxide emission in the J = 2 — 1 
transition of the three brighter isotop ologues (^^CQ/ '^CQ,C^^Q) towards NGC6000 



using the Submilhmeter Array (SMA; 



Hoetal 



2004J ). This is the first high resolution 



morphological study of the molecular component towards NGC 6000. 



2. Observations 

Observations were carried out on 2007 February 20 with the Submilhmeter Array 
(SMA) atop Mauna Kea, Hawaii. Seven out of the eight antennas of the array were 
operational and placed in the "compact-north" configuration, with baselines ranging 
16 — 120 m. The SIS receivers were tuned to the CO J = 2 — 1 frequency (230.538 GHz) 
in the upper side band (USB), re dshifted to the NGC 6000 average systemic velocity of 



2145 km s"^ (lAlbrecht et al 



20071 ). The correlator was configured so that the J = 2 — 1 
transitions of the isotopic substitutions ^^CO (220.398 GHz) and C^^O (219.560 GHz) were 
simultaneously observed in the lower side band (LSB). Spectral resolution was 0.8125 MHz 
(~ Ikms^^). System temperatures ranged from 125 — 250 K, with a zenith opacity of 
0.1 — 0.15 as measured at 225 GHz from CSO. 

A single field (field of view ~ 55" at 230 GHz) was observed towards th e central positi on 
g j^nn n n = 15''49'"49!40 and 5j2ooo.o = -29°23'11'.'0 as measured by IRAS Jsanders et al.l 



20031 ) . We achieved a resolution of 2.8" x 1.7" in the uniform weighted maps and 3.7" x 2.2" 
in the natural weighted maps. At the distance of NGC 6000, 1" is equivalent to a projected 
distance of ~ 153 pc. Passband calibration was derived from spectra of Saturn and 3C 279. 
Absolute flux calibration (within a 10% accuracy) was determined from observations of 
Callisto. Gain fluctuations were corrected with the nearby quasar J1517-2422 at a distance 
of 8.7° from NGC 6000. The applied gain corrections were checked by imaging another 
nearby quasar, J1625-2527, at a distance of 8.9°. Data calibration and reduction was 
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performed using the MIR-IDL package and imaging using MIRIAD (jSault et al.lll995l ). 



Results 



3.1. Continuum Emission 



Continuum emission from the nucleus of NGC 6000 was imaged using the CO emission 
free channels in the USB. The spectral energy distribution (SED) of NGC 6000 shows 
how most of its continuum emission at 1.3 mm stems from thermal dust emission with 



no significant evidence of free-free contribution ( iRoche fc Chandler 



19931). In Fig. m {left 



panel) we show the 1 . 3 mm continuurn emis sion over the an average image of the J, H 



20061 ). The position of the peak of emission at 



and K 2MASS bands (ISkrutskie et al. 
aj2ooo.o = 15^49°^49'!5, 5j2ooo.o = -29°23'13" matches that of 2MASS emission within the 
resoluti on. This position als o matches within 2" that of the FIR from low resolution IRAS 



images (jSanders et al. 



Condon et al. 



2003|) and within 1" that of the VIA 1.4 GHz peak (18" resolution, 



19961 ). Together with the CO mapping presented in Sect. 13.21 we have 



obtained the most accurate determination of the nucleus position in NGC 6000, with a 



resolution better than the 2.5" of 2MASS (ISkrutskie et al. 



20061 ). 



The faint unresolved continuum source has a total flux of 13 ± 3 mJy as measured 



from our observations. This flux seems to agree wit h the 18 ± 8 mJy f 



1.3 mm with the 19.5" beam of the JCMT telescope (IRoche &: Chandler 



ux observed at 



19931 ). However, 



the JCMT flux might be overstimated given that the contribution of the CO line (Sect. 13. 2p 
to the continuum in the 64 GHz bolometer band would be ~ 13 mJy. Therefore, the 
corrected JCMT continuum flux would be a factor of two lower than our measurement. 
On the other hand, the CO emission would contribute ~ 16mJy to the 42.7 ± 3.9 mJy 
average flux measured by the 11" beam of the IRAM 30 m over its 50 GHz bolometer band 
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( IChini et al.lll995l ). This corrected flux is a factor of two larger than what we recover from 
our observations. An even larger 1.3 mm flux of 60.6 ± 7.8 mJy was measured by SEST 
telescope, which is likely d ue to its larger beam (24") collecting most of the extended 



emission ( Chini et al. 



19951 ). Indeed, the extended 1.3 mm continuum emission is estimated 



to be 98mJy (CO c orrected) in the inner 70" with the 5 pointing map observations with 



the SEST telescope flAlbrecht et al. 



20071 ). Assuming the continuum measurements from 



SEST, we only recovered the compact ~ 15% of the continuuum emission in our maps. 



3.2. Molecular Gas Emission 

3.2.1. Comparison to single dish 
The total integrated flux recovered from the ^^CO 2 — 1 image is 1080 ± 60 Jykms"^. 



Single dish observations of this transition with the 24" beam of t 



in a total detected flux of 780 ± 80 Jy km s ^ (iChini et al. 



le SEST telescope resulted 



19961 ). where we assumed a 



conversion factor of ~ 25 Jy/K. This difference is contrary to the expected lower flux in 
the interferometric maps due to filtering of extended emission. Indeed, the ~ 30% flux 
difference is too large to be attributed to absolute flux calibration error of both SEST and 
SMA data. Although SEST observations were taken towards a nominal position 2" away 
from the position of peak emission derived from our map, this difference would only explain 
a < 4% lower measured flux. In order to understand this measured fl ux, we convo l ved o ur 



Chini et al. 



f ll996h 



map to a 24" resolution. The asymmetric lineshape in the data from 
clearly differs from that obtained at the center of the convolved map and significantly 
matches the lineshape at a position > 5" south of the central position, consistent with the 
SEST pointing accuracy. This difference in the observing position easily accounts for a 
difference in the integrated intensity of > 20% with respect to the central position. Though 
the uncertain comparison with single-dish data does not allow us to accurately determine 
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the missed flux in our maps, we can assume that no significant amount of flux have been 
filtered out and that most of the extended CO emission has been recovered in our maps. 



3.2.2. CO morphology 
The total integrated ^^CO J = 2 — 1 emission detected in the nuclear region of 



NGC 6000 is presented in Fig. [T] {right panel) on top of the optical HST image ( ICaroUo et al. 



19971 ). In addition to the strong nuclear emission, a fainter extended structure is clearly 
detected (> 6a) extended in the north-south direction, up to a distance of ~ 12" (1.8 kpc) 
from the center. A similar north-south structure is observed in the near-infrared (Fig. [H 

left panel). 

In Fig. [2] we present the maps of the CO J = 2 — 1 integrated emission over the range 
from 1900 to 2500 km s~^ {Left panel) and the intensity weighted mean velocities {Right 
panel). The bright nuclear region is just barely resolved by the 2.8" x 1.7" (430 x 260 pc) 
beam. A Gaussian fit to the UV visibilities in the nuclear emission results in a deconvolved 
extent of 3.7" x 2.1" (560 x 320 pc), centered at ajaooo.o = 15M9'"49!5, 5j2ooo.o = -29°23'13" 
(see Table E]). If we assume that the deconvolved ellipticity is due to the inclination of the 
circumnuclear molecular disk (Sect. H]), an inclination of 31° is inferred. This value is in 
agreement with those derived from large scale images of NGC 6000. The observed position 
of the CO emission peak matches the derived dynamical center and is in agreement with 
that derived from the continuum peak. Similar to what is observed in the near-IR 2MASS 
images, the peak of CO emission does not coincide with the optical peak of emission due to 
the strong extinction towards the nuclear region. In fact, the extended CO emission follows 
the dust lanes seen in the optical. 

Optical imaging of the nuclear region of NGC6000 shows a bar-like structure with a 
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nucleus surrounded by aii irregular star forming ring with spiral arms down to the ring 



( ICaroUo et al. 



1997 



20021 ). Fig. [3] shows the red and blue shifted emission in a close-up 
view of the central region. The star forming ring structure, offset by {a, 6) ~ ( — 1",!") from 
the peak of molecular emission, appears to be located at the peak of the high velocity CO 
emission (Fig. [3]), right at the eastern edge of the molecular disk. This structure is probably 
located closer to the observer than the gas disk, and therefore less affected by obscuration. 
The strongly asymmetrical optical image with respect to the galactic center points to a 
high extinction, as derived in Section I3.3.2[ likely hiding the nuclear and eastern region. 



3.2.3. CO kinematics 

Fig.HJshows the integrated CO emission over 30kms~^ channels. The resolved molecular 
emission in the nuclear region is observed over a range of velocities Av > 400 km s^^ while 
the more extended emission is concentrated in the central Av ~ 300 km s^^. As shown in 
the channel maps (Fig. H]) the optical bright structures are located at the regions where the 
most red shifted emission {v > 2300 km s~^) peaks. 

The mean velocity map in Fig. [2] [Right panel) shows a strong velocity gradient in 
the central region, with a smoother variation outside the nucleus. This velocity gradient 
is clearly shown by the position-velocity (P-V) diagrams shown in Fig. [51 P-V diagrams 
have been plotted in two different directions, first across the whole molecular emission 
in the north-south direction (P.A.=0°; left panel in Fig. |5]) and second in a direction 
perpendicular to the axis of rotation of the nuclear region (P.A.=133°; right panel in Fig. |5]). 
Both directions are indicated in the right panel in Fig. [2] as continuous and dashed lines, 
respectively. The P-V diagrams in both directions show how in the central 4" (600 pc) the 
molecular gas disk is rotating like a rigid-rotor with a rotation velocity ~ llOkms"^ per 
arcsec, equivalent to 0.7 km s^^ pc^^, after correcting the velocity for an inclination of ~ 30° 
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f lTullv 



1988 



Corwin et al 



the sample of galaxies by 



19851). This rotation velocity is similar to the average found in 



Sakamoto et al. 



(Il999l ) containing barred galaxies like NGC 6000. 



3.2.4. CO and C^^O emission 



18/ 



The emission of the rarer isotopologues were detected towards the nuclear region. 
Fig. [6] shows the integrated emission of ^'^CO and C^^O compared to that of the main 
isotopologue. These maps have been natural weighted in order to get a higher sensitivity at 
the expense of a larger 3.7" x 2.2" (560 x 340 pc) beam. The nuclear emission is resolved in 
both ^^CO and C^^O in the direction of the disk rotation. However, no extended emission 
in the north-south direction has been detected in the rarer species due to their much lower 
abundance and lack of enough sensitivity. Still, the two lowest contours plotted for ^^CO 
appear to be slightly elongated in the same direction as the extended emission traced by 
^^CO. Gaussian models fitted to the visibilities of each isotopologue give the parameters 
listed in Table [2l We observe the emission of all the three isotopical substitutions to peak 
at the same position. The spectrum at the central synthesized beam of each of the natural 
weighted datacubes is shown in Fig. U\ where the ^'^CO and C^^O spectral features have 
been multiplied by a factor of 10 for the sake of comparison with the much brighter CO 
feature. 



From our observations we derive an integrated intensity ratio -R2-1 ~ -^[^^CO j=2-i]/-^[^'^C0 



J=2-lJ 



of ~ 15.5. This ratio is in agreement with those Rl 



13 



11.3 ±3.3 and R^i^ 



12.7 ±4.7 



ratios derived from the CO J = 1 — and J = 2 — 1 data compilations on n ormal 
galaxies with Lpm < 10^^ Lq ( iTaniguchi fc Ohyamal . 



1998 



Taniguchi et al 



19991 ). The ratio 



-R2-1 ~ 57.9 calcul ated with the C O isotopolog ue is similar to those derived for other 



19931 ). Therefore, assuming the isotopic ratios 



starburst galaxies (iHenkel &: Mauersbergerl 
of i2c/i3c= 40 - 50 and ^^0/^^0= 150 - 200 apply to NGC 6000, we can estimate an 



10 



opacity of the main isotopologue line of tco{2-i) ~ 2.5 — 5. We can then accurately assume 
both ^^CO and C^^O to be optically thin. 



With the derived parameters in Table [2] for central spectra of the CO map we can make 
a rough estimate of the total dynamical mass within the nuclear region of NGC 6000. The 
dynamical mass can be calculated as a function of the disk size and the maximum velocity 
corrected by inclination as Mdyn = 250 Mq {Vmax/ sin{i)y R^, where the Vmax and -Rd are 
given in kms~^ and pc, respectively. For this calculation we used the peak velocity of 
Vmax ~ 180kms~^, with respect to the systemic velocity, reached at a distance of Rd ~ 1.5" 
as measured in the south-east direction of the molecular disk (low velocity end in Fig. |5]) so 
that we avoid the uncertainty due to the disk asymmetry at the higher velocities (Sect. SJ. 
This approximation yields a total mass M^yn = 7 x 10^ Mq. 



In order to get an accurate estimate of the mass of molecular hydrogen traced by the 
carbon monoxide we will use the total flux detected of the isotopologue C'^^O emission. 
This is equivalent to using the main isotopologue corrected by opacity effects. Assuming 
optically thin emission and LTE conditions we can get the mass of hydrogen as 



where ^^C^^O/^C^O is the ratio between the main form of carbon monoxide and any 
given isotopologue. This expression can be applied for the three isotopologues and 



3.3. 



Mass determination 



3.3.1. Dynamical mass 



3.3.2. Mass of molecular gas traced by CO 




(1) 
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assumes a CO fractional abundance ^^C0/H2 ~ 8.5 x 10 ^ 



Tex = 30 ± 10 K, similar to the dust temperature derived by 



(IFrerking et al. 



Yang fc Phillips 



Paglione et al 



982) 



and 



(120071). Under 



{2001)^With the total 



1999I ). we derive an 



these conditions, Eq. [T] is equivalent to Eq. 3 in 
flux observed in C^^O and the ratio ^^0/^®0~ 150 (IHarrison et al. 
average A^H2 = (1-0 ± 0.3) x lO^^cm"^. This amount of molecular hydrogen implies a 
visual extinction of ~ 100 magnitudes towards the central region of NGC 6000 which is 
certainly responsible for the asymmetric optical structure of the nucleus of NGC 6000. 

We estimate a total mass of molecular gas in the nuclear region of Mh2 ~ 
2.9 ± 0.9 X 10® Mq within the central projected ~ 566 x 321 pc (as measured from the 
deconvolved CO source size). 

From the CO map in Fig. [T] we estimate that ~ 85% of the molecular mass in NGC 6000 
is concentrated in the nuclear ~ 3" region (~ 460 pc). This percentage is just a lower limit 
given the significant opacity affecting the CO emission mostly towards the central peak of 
emission. Therefore we can estimate a total mass of molecular gas of 3.5 x 10® Mq within 
the inner ~ 1.8 kpc. Assuming a correction factor 1.36 for the He and heavier elements 



(Allen 



19731 ) ■ we can estimate a total mass of gas Mgas = 4.6 x 10® Mq, which yields a ratio 



Mgas/Mayn ~ 6%. 



4. Discussion 

4.1. Bar-driven molecular gas fueling the starburst in NGC 6000 



The study of the molecular gas content in a sample of 20 galaxies by 



Sakamoto et al. 



(119991 ) statistically shows that the gas tends to be more concentrated in the central 
kiloparsec in barred systems. This is the case of NGC 6000 where most of the gas is located 
in the inner half kiloparsec (Sect. 13.3.21) . The large-scale bar revealed in the NIR may also 
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have a fingerprint in the molecular gas kinematics described in this work. We observe that 
the P.A.=0° P-V diagram in Fig. O (right) shows how the extended molecular component 



displays a S-shaped feature. This kinematic s ignature can 



De understood as the noncircular 



motions in the context of a barred potential (IBinney et al. 



199ll ). The inner material could 



be moving in the barely resolved X2 orbits which, with the resolution of our maps, would 
look like a molecular disk. On the other hand the external gas would be moving in large 
elliptical xi orbits traced as the S-shaped profile in the P-V diagram. Higher resolution 
imaging of NGC 6000 might support the scenario of a barred potential b y resolving the inner 



disk structure into circular X2 orbits as observed by 



Meier et al 



(120081 ) towards Maffei2. 



Similar signatures of a barred pot ential are observed in the P-V diagrams o f other galaxies 



such as NGC 1530 and NGC 4258 f Pownes et al. 



1996 



Cox fc Downes 



19961 ) and equivalent 



structure are found in their innermost regions. 

However our mea sured Mgas/Mdyn ratio of 6% is lower than the average measured by 



Sakamoto et al. 



( 119991 ) in starbursts and barred galaxies. This measurement is significantly 



affected by the way the gas mas is calculated Indeed, adopting their same conversion 



factor (X = 3 X 10^^ cm ^(Kkms ^) ^ rather than X = 0.4 x 10^^ cm ^(Kkms 



see 



Sect. 14. 3p the measured ratio would increase u p to ~ 35%. This ratio is only found in 



barred starbursts in the 



Sakamoto et al. 



fll999h sample. 



4.2. The high velocity asymmetry 

The strong asymmetry in the nuclear ring/disk is also a particularly interesting feature. 
As seen in Fig. [5] a significant part of the molecular gas traced by CO is skewed towards the 
region at the lowest velocities, corresponding to the south-east disk component. However, 
the emission at the highest velocities from 2300 to 2400 km s~^ in the P-V diagrams, 
although following the same rotation gradient, seems to be significantly detached from the 
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inner rotating structure. This feature at the highest velocities corresponds to the wing 
observed in the spectrum shown in Fig. [71 We discard the possibihty of self- absopt ion at 
high velocities, which would change the systemic velocity derived from the Gaussian fit to 
the line, as the optically thin ^^CO and C^^O show similar profiles and are centered at 
the same velocity. It is surprising that ^'^CO and even are detected in this velocity 

component between 2300 to 2400 km s~^. Moreover, the low ratio of ^^CO and C^^O 
with respect to the main isotopologue implies that a significant opacity is affecting even 
the ^^CO emission in this molecular component. CO being optically thick can provide 
some constraint on the extent of the emission as compared to our resolution. From the 
CO brightness temperature measured Tb ~ 1 K and assuming an excitati on temperature 
Tex = 30 K (as derived from the dust temperature, lYang fc Phillips! 120071 . and assuming 
CO being thermalized at this temperature), we estimate the region extent to be < 80 pc 
(< 0.5"). This is slightly above the ~ 50 pc typical size of giant molecular clouds within 
our Galaxy. This giant molecular cloud within the inner region of NGC 6000 must be very 
dense and more compact than our estimate to explain such CO opacity. It could have 
originated as a shock between the circumnuclear disk and the infalling molecular gas along 
the dust lanes. Moreover, emission arises from the same region as the star forming ring 
which suggest that this dense and compact molecular cloud might be directly related to the 
fueling of star formation event in this region. 



4.3. CO-to-H2 conversion factor in starburst galaxies 



We can estimate the column density of molecular gas via the conversion factor 
X = N{H2) / CO = 3. X 10^° cm'^fK km s~^)~^ based on measurements of galactic 



disk molecular clouds (ISolomon et al 



19871 ). This way, and adopting the CO integrated 



intensity in this work, we derive a H2 column density of 7.5 x lO^^cm ^ which is almost 
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an order of magnitude larger than the value derived from the optic ally thin C^^O line 



(Sect . I3.3.2p. As already o bserved towards the starburst NGC 253 flMauersberger et al. 



199 



Harrison et al 



19991 ) ■ the conversion factor in the starburst environment might be 
significantly lower than that in the Galactic disk. From our measurements we calculate 
a conversi on factor X = 0.4 x 10^^ cm~^(Kkms~^)~^ in agreement with that derived for 



NGC 253 flMauersberger et al. 



19961 ). This result support the evidences found for a lower 



conversion factor in t he central region of galaxi es, and in particular in starburst with respect 



to the Galactic disk (jPownes fc Solomon 



the Large Magellanic Clouds ( IWang et al. 



9981) or that measured in star forming regions in 



20091). 



5. Summary 

We have obtained high resolution (2.8" x 1.7") maps of the molecular emission 
towards NGC 6000 as traced by the J = 2 — 1 transition of carbon monoxide with the 
SMA. Emission from the three main isotopologues of carbon monoxide, ^^CO, ^^CO, 
and C^^O are detected. Both the continuum and molecular emission peak at the same 
position. From these measurements we have accurately determined the nucleus position in 
NGC 6000 at aj2ooo.o = 15'^49™4g:5 and 5j2ooo.o = -29°23'13, in agreement with previous 
IR determinations. More than an 85% of the emission is concentrated towards the nuclear 
projected ~ 400 x 300 pc region, where we estimate optical extinctions of the order of ~ 100 
magnitudes. Such large obscuration is likely responsible for the asymmetry in the optical 
images, probably hiding the eastern part of the nuclear region at these wavelengths. 

The emission from the central region is resolved in the direction of the observed 
velocity gradient. Molecular emission is detected in the north-south direction extended 
up to distances of > 1.5 kpc from the center. However this is only detected in the main 
isotopologue. Both the nuclear and extended emission show dynamical evidences that might 
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be explained by the presence of a bar potential affecting the molecular material, likely 
responsible of the strong nuclear concentration and maybe the starburst event. We observed 
a high velocity component which appears to be detached from the nuclear disk/ring and is 
highly optically thick. With an estimated extent of < 80 pc, this giant molecular cloud is 
likely to be related to the fueling of the nuclear star burst in NGC 6000. 

The ratios of ^^CO and C^^O with respect to the main isotopologue are similar to the 
average m easured in similar starburst g alaxies with luminosities Lpm < lO^^L© in the 



literature (iHenkel fc Mauersberger 



19931) 



The total dynamical mass derived for NGC 6000 is 1.9 x 10^ Mq. We estimate a total 
mass of molecular gas of 4.6 x lO^ Mr;^ and a ratio Mgns/M dyn ~ 6%, similar to the starbursts 



in the barred galaxy sample by (ISakamoto et al. 



19991 ). The independent estimate 



of the molecular mass with the optically thin C^^O isotope allows the determination 
of the conversion factor C0-to-H2 of Xco = 0.4 x 10^^ cm~^(Kkms~^)~^, which 
agrees with previous estim ates that measure a lower factor in starburst nuclear regions 



(IMauersberger et al. 



19961 ) as compared with the standard factor derived for Galactic disk 



clouds. 



This paper uses data taken with the Submillimeter Array on Director's discretionary 
time as part of the Astronomy 191 course for Harvard University undergraduates. We 
want to thank Prof. Carl Heiles for his valuable comments on the manuscript. The 
Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory 
and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the 
Smithsonian Institution and the Academia Sinica. 

Facilities: SMA () 
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contours SMA CO 2-1 
grey HST F606W 




IS'Mg^SO^'O 49=^5 49''0 
RA (J2000) 



IS'Mg^aO^'O 49=^5 49''0 
RA (J2000) 



Fig. 1. — Left panel: 1.3 mm continuum natural weighted emission in NGC6000 as 



measured with the SMA 
scale 



Skrutskie et al. 



'contours) over an average of the J, H, and K 2MASS (grey 
20061 ). Contours levels are 2a significant as —5, 5, and lOmJy. 
2MASS image is displayed in logarithmic scale. Right panel: Composite image of the 
uniform weig hted 12 CO J = 2 - 1 integrated emi ssion (contours) ari d the optical HST- 



WFPC2 image in the F606W filter of NGC6000 flCaroUo et al.lll997f ). Contours are 3a 
levels from 12 to 74 Jy beam"^ kms"^ and lOcr from 115 to 445 Jybeam"^ kms"^, where 
a = 4.1Jy/beamkms~i in the image. HST image is displayed in logarithmic scale. The 
respective synthesized beams of each map are displayed in the bottom left-hand corner of 
each plot. 
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15^49"'50"0 49"5 49!'0 50'0 49^5 49^0 

^J2000 ^J2000 



Fig. 2. — ^^CO J = 2 ^ 1 integrated emission and kinematics towards the nuclear region of 
NGC6000. The synthesized beam of the uniform weighted maps (2.8" x 1.7") is plotted at 
the bottom left-hand corner of each plot. Left panel: Integrated flux over the velocity range 
of 1900 and 2500 km s~^. Contours are similar to those used in Fig. [TJ Right panel: Mean 
velocity field map. Contours are uniformly spaced by 50kms~^. Only pixel above 3a level in 
integrated intensity map are shown. The solid and dashed lines are the cut directions used 
for the position-velocity diagrams in Fig. [51 
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Fig. 3. — Close-up image of the central region of NGC6000. The CO red (solid contours) 
and blue shifted (da shed contours) emission are plotted over the optical HST-WFPC2 image 



(jCarollo et al. 



19971 ). Emission has been integrated —250 and +250 km s ^ with respect to 



the systemic velocity (2135 km s~^) for the blue and red shifted components, respectively. 
Contours are 3(j levels from 7.5 to 60Jykms~^ and lOcr from 60 to 185Jykms~^, where 
a ~ 2.5 Jy/beamkms~^ in both images. The optical emission is strongly asymmetrical 
with respect to the dynamical center. We observe that the star forming ring, significantly 
offset from the center of the galaxy, is located at the position where the highest velocity CO 
emission is observed as better seen in the channel maps in Fig. HI 
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Fig. 4.- 



12 



CO channel map s (Contours) in steps of 30kms ^ plotted over the HST-WFPC2 



image ( ICaroUo et al. 



19971 ). Contours are plotted at 2a levels of l.SJykms ^. Central 



velocity of each channel is shown in the upper left corner and beam size in the lower right. 
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Fig. 5. — Position- Velocity diagrams across the nuclear region of NGC6000 in the P.A.=0° 
and P.A. = 133° directions. Contours are 6a levels (0.15 Jy) from 0.15 Jy to 1.95 Jy. Both 
cuts directions are shown in Fig. [2] by the solid and dashed lines, respectively. 
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Fig. 6. — Natural weighted integrated flux maps of the observed J = 2 — > 1 transition in 
the different CO isotopologues (contours), compared to the ^^CO emission (grey scale). The 
beam size in these maps is 3.7" x 2.2" (shown as filled ellipses at the bottom left corner). ^^CO 
grey scale levels corresponds to those in Fig. [2j Contours are 2a levels (4.4 Jybeam"^ kms^^) 
in both i^CO and C^^O. 
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Fig. 7. — Brightness temperature spectra of the three isotopologues extracted from the 
central synthesized beam of the natural weighted datacubes. The resolution in velocity is of 
lOkms"^ for the main isotope, 20kms~^ for ^^CO, and 30kms~^ for C^^O. The spectra of 
^^CO and C^^O has been multiphed by a factor of 10 for the sake of clarity in the comparison. 
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Table 1. Summary of observational parameters for NGC 6000 



Parameter Value Reference 



GALAXY POSITION 

ROSAT (J2000) 15:49:50.1, -29:23:10.7 1 

IRAS (J2000) 15:49:49.4, -29:23:11 2 

Flsr 2145kms-i 3 

Distance (iJo = 75 A:TOs-iMpc-i) 31.6Mpc 4 

MORPHOLOGY 

Morphological type SB(s)bc 5 

Inclination 33°, 27° 6,7 

logD25(0.1') 1.27 5 

Weighted mean luminosity class 4.0 7 

T^ldTTiGt-GfjyiQ^ j jnin 2.3x2.0 8 

RA. 154°, 163° 5,9 

PHYSICAL PARAMETERS 

Trf„,t(K) 34,30 10,11 

Mdust {Mq) 1.5 X 10^, 2.4 X 10^ 10,11 

log Lfir{Lq) 10.83 12 

log LiRiLQ) 10.97 12 

log Mhi (Mq) 9.86, 9.8 13,14 

log Mh2 {Mq) 8.5 15 

log Mgas {Mq) 8.7, 10.0 15,16 

log Mdyn {Mq) 9.8 15 
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Table 2. Derived parameters from the CO, ^^CO and C^^O observed emission 



Paraniotcr CO ^'^CO C^^O 

Gaussian fit to the maps ^ 

Center (J2000.0) 15:49:49.5, -29:23:13 15:49:49.5, -29:23:13 15:49:49.5, -29:23:13 

Size (FWHM) 3.7" x 2.1", P.A. 12° 2.7" x 2.1", P.A. 29° 2.5" x 1.3", P.A. 5° 

Total Flux 922 ±9 Jy km s"^ 90 ± 7 Jy km s"^ 33 ± 6 Jy km s"^ 

Central pixel spectra ^ 

Vlsr 2135.3 ± 0.3 kms"^ 2119±4kms-i 2104±9kms-i 

At;i/2 217.1 ± 0.6 kms-i 196±10kms-^ 150±20kms-i 

Tb 5600 mK 400 mK 140 mK 



^Results of the deconvolved Gaussian fit models applied to the UV visibihties. 
^Parameters from the Gaussian fit to the spectra at the central pixel of each map. 



-25 - 

REFERENCES 

Albrecht, M., Kriigel, E., & Chini, R. 2007, A&A, 462, 575 
Allen, C. W. 1973, Astrophysical quantities, ed. C. W. Allen 

Binney, J., Gerhard, O. E., Stark, A. A., Bally, J., & Uchida, K. 1. 1991, MNRAS, 252, 210 

Boiler, T., Meurs, E. J. A., Brinkmann, W., Fink, H., Zimmermann, U., & Adorf, H.-M. 
1992, A&A, 261, 57 

CaroUo, C. M. 1999, ApJ, 523, 566 

CaroUo, C. M., Stiavelli, M., de Zeeuw, R T., & Mack, J. 1997, AJ, 114, 2366 

CaroUo, C. M., Stiavelli, M., Seigar, M., de Zeeuw, P. T., & Dejonghe, H. 2002, AJ, 123, 

159 

Chini, R., Kruegel, E., & Lemke, R. 1996, A&AS, 118, 47 

Chini, R., Kruegel, E., Lemke, R., & Ward-Thompson, D. 1995, A&A, 295, 317 

Condon, J. J., Helou, C, Sanders, D. B., & Soifer, B. T. 1996, ApJS, 103, 81 

Corwin, H. C, de Vaucouleurs, A., & de Vaucouleurs, G. 1985, University of Texas 
Monographs in Astronomy, 4, 1 

Cox, P., & Downes, D. 1996, ApJ, 473, 219 

de Vaucouleurs, C, de Vaucouleurs, A., Corwin, Jr., H. C, Buta, R. J., Paturel, C, & 
Fouque, P. 1991, Third Reference Catalogue of Bright Galaxies (Volume 1-3, XII, 
2069 pp. 7 figs.. Springer- Verlag Berlin Heidelberg New York) 

Downes, D., Reynaud, D., Solomon, P. M., & Radford, S. J. E. 1996, ApJ, 461, 186 



-26- 



Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 

Dressier, A. 1991, ApJS, 75, 241 

Fathi, K., & Peletier, R. F. 2003, A&A, 407, 61 

Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590 

Harrison, A., Henkel, C, & Russell, A. 1999, MNRAS, 303, 157 

Henkel, C, & Mauersberger, R. 1993, A&A, 274, 730 

Ho, R T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 

Hunt, L. K., & Malkan, M. A. 2004, ApJ, 616, 707 

Jogee, S., Scoville, N., & Kenney, J. D. P. 2005, ApJ, 630, 837 

Knapen, J. H., Perez-Ramirez, D., & Laine, S. 2002, MNRAS, 337, 808 

Koribalski, B. S., Staveley-Smitli, L., Kilborn, V. A., Ryder, S. D., Kraan-Kortcwcg, R. C, 
Ryan-Weber, E. V., Ekers, R. D., Jerjen, H., Henning, P. A., Putman, M. E., Zwaan, 
M. A., de Blok, W. J. G., Calabretta, M. R., Disney, M. J., Minchin, R. F., Bhathal, 
R., Boyce, P. J., Drinkwater, M. J., Freeman, K. C, Gibson, B. K., Green, A. J., 
Haynes, R. F., Juraszck, S., Kcstcven, M. J., Knczck, P. M., Mader, S., Marquarding, 
M., Meyer, M., Mould, J. R., Oosterloo, T., O'Brien, J., Price, R. M., Sadler, E. M., 
Schroder, A., Stewart, I. M., Stootman, F., Waugh, M., Warren, B. E., Webster, 
R. L., & Wright, A. E. 2004, AJ, 128, 16 

Martin, J. M., BottineUi, L., Gouguenheim, L., & Dennefeld, M. 1991, A&A, 245, 393 

Mauersberger, R., Henkel, C., Walsh, W., & Schulz, A. 1999, A&A, 341, 256 



-27- 

Mauersberger, R., Henkel, C, Wielebinski, R., Wiklind, T., & Reuter, H.-R 1996, A&A, 
305, 421 

Meier, D. S., Turner, J. L., & Hurt, R. L. 2008, ApJ, 675, 281 

Paglione, T. A. D., Wall, W. F., Young, J. S., Heyer, M. H., Richard, M., Goldstein, M., 
Kaufman, Z., Nantais, J., & Perry, G. 2001, ApJS, 135, 183 

Perez- Ramirez, D., Knapen, J. H., Peletier, R. F., Laine, S., Doyon, R., & Nadeau, D. 2000, 
MNRAS, 317, 234 

Pfenniger, D., & Norman, C. 1990, ApJ, 363, 391 

Pizzella, A., Corsini, E. M., Dalla Bonta, E., Sarzi, M., Coccato, L., & Bertola, F. 2005, 
ApJ, 631, 785 

Roche, P. F., & Chandler, C. J. 1993, MNRAS, 265, 486 

Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999, ApJ, 525, 691 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C., Surace, J. A., & Soifer, B. T. 2003, AJ, 126, 
1607 

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in ASP Conf. Ser. 77: Astronomical 
Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, & J. J. E. 
Hayes, 433 

Sheth, K., Vogel, S. N., Regan, M. W., Thornley, M. D., & Teuben, P. J. 2005, ApJ, 632, 
217 

Siebenmorgen, R., Kriigel, E., & Spoon, H. W. W. 2004, A&A, 414, 123 



-28- 



Skrutskie, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., Schneider, S., Carpenter, 
J. M., Beichman, C, Capps, R., Chester, T., Ehas, J., Huchra, J., Liebert, J., 
Lonsdale, C, Monet, D. G., Price, S., Seitzer, P., Jarrett, T., Kirkpatrick, J. D., 
Gizis, J. E., Howard, E., Evans, T., Fowler, J., Fullmer, L., Hurt, R., Light, R., 
Kopan, E. L., Marsh, K. A., McCallon, H. L., Tarn, R., Van Dyk, S., & Wheelock, 
S. 2006, AJ, 131, 1163 

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 

Taniguchi, Y., & Ohyama, Y. 1998, ApJ, 507, L121 

Taniguchi, Y., Ohyama, Y., & Sanders, D. B. 1999, ApJ, 522, 214 

TuUy, R. B. 1988, Nearby galaxies catalog (Cambridge and New York, Cambridge University 
Press, 1988, 221 p.) 

Wang, M., Chin, Y.-N., Henkel, C, Whiteoak, J. B., & Cunningham, M. 2009, ApJ, 690, 
580 

Yang, M., & Phillips, T. 2007, ApJ, 662, 284 

Young, J. S., Xie, S., Tacconi, L., Knezek, P., Viscuso, P., Tacconi-Garman, L., Scoville, 
N., Schneider, S., Schloerb, F. P., Lord, S., Lesser, A., Kenney, J., Huang, Y.-L., 
Devereux, N., Claussen, M., Case, J., Carpenter, J., Berry, M., & Allen, L. 1995, 
ApJS, 98, 219 



This manuscript was prepared with the A AS macros v5.2. 



